Searching for the fastest dynamo: 
Laminar ABC flows 
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The growth rate of the dynamo instability as a function of the magnetic Reynolds number is 
investigated by means of numerical simulations for the family of the ABC flows and for 2 different 
forcing scales. For the ABC flows that are driven at the largest available length scale it is found 
that as the magnetic Reynolds number is increased: (a) The flow that results flrst in dynamo is the 
2iD flow for which A=B and C=0 (and all permutations), (b) The second type of flow that results 
in dynamo is the one for which A = i? ~ 2C/5 (and permutations), (c) The most symmetric flow 
A=B=C is the third type of flow that results in dynamo, (d) As Rj^^ is increased, the A=B=C flow 
stops being a dynamo and transitions from a local maximum to a third-order saddle point, (e) At 
larger ii^^ the A=B=C flow re-establishes its self as a dynamo but remains a saddle point, (f) At 
the largest examined the growth rate of the 2iD flows starts to decay, the A=B=C flow comes 
close to a local maximum again and the flow A = B ~ 2C/5 (and permutations) results in the 
fastest dynamo with growth rate 7 ~ 0.12 at the largest examined Rj^^ . For the ABC flows that are 
driven at the second largest available length scale it is found that (a) the 2^0 flows A=B, C=0 (and 
permutations) are again the flrst flows that result in dynamo with a decreased onset, (b) The most 
symmetric flow A=B=C is the second type of flow that results in dynamo. It is and remains a local 
maximum, (c) At larger Rm the flow A = B ~ 2C/5 (and permutations) appears as the third type 
of flow that results in dynamo. As 7?^^ is increased it becomes the flow with the largest growth rate. 
The growth rates appear to have some correlation with the Lyaponov exponents but constructive 
re-folding of the field lines appears equally important in determining the fasted dynamo flow. 



I. INTRODUCTION 



Magnetic dynamo is the process through which an elec- 
trically conducting fluid amplifles and maintains mag- 
netic energy against Ohmic dissipation by continuously 
stretching and re-folding the magnetic field lines [l|, 
This process is considered to be the main mechanism for 
the generation of magnetic energy in the universe. It 
is present in the intergalactic and interstellar medium, 
in accretion disks and in the interiors of stars and plan- 
ets. It has also been realized recently in different labora- 
tory experiments 0-111 • The flows in these examples vary 
in structure and the generated magnetic fields exhibit a 
large variety of structural and temporal behavior. It is 
then desirable to understand which properties of a flow 
are important for the amplification of magnetic energy 
and how do they effect the dynamical behavior of the 
magnetic field. This question is of particular interest for 
the dynamo experiments for which optimizing the flow is 
important for achieving dynamo at small energy injection 
rates [1]. 

In theoretical studies, various flows have been exam- 
ined analytically and numerically both in the laminar and 
in the turbulent regime. The Ponomarenko 0, the ABC 
the Roberts fnl, [13], the Taylor-Green [11], and 
the Archontis flow [ij] are some of the flows that have 
been shown to result successfully in dynamo action pro- 
vided that the magnetic Reynolds number R^^j (the ratio 
of the large-scale velocity time scale to the large-scale dif- 
fusivity time scale) is sufficiently large. The choice of flow 
for study was motivated either by its similarity to astro- 



physical flows or due to its simplicity that allowed analyt- 
ical treatment or made the investigation more tractable 
numerically. Other than this practical motivation there 
is no mathematical justiflcation for preferring one flow 
over an other. 

This lack of mathematical reasoning motivates this 
work. Over a family of flows of flnitc energy and vortic- 
ity not all members are as efficient in producing dynamo 
action. It is then expected that a fiow in this family 
exist that is optimal for dynamo action. Finding and 
investigating the properties of such an optimal fiow can 
then reveal which mechanisms are important for mag- 
netic field amplification. How an optimal fiow is defined 
is described in the next section where the general prob- 
lem is formulated in detail. 



II. FORMULATION 

At the early stages of the dynamo, when the Lorentz 
force is too weak to act back on the fiow, the evolution 
of the magnetic field is given by the linear advection dif- 
fusion equation 



9tb + u • Vb = b • Vu + T^V^b, 



(1) 



where b is the magnetic vector field, u is the velocity field 
and 7] is the magnetic diffusivity. The advection term in 
the left hand side of Eq. ([T|) is responsible for the mixing 
of the magnetic field lines. The first term in the right 
hand side is the stretching term that is responsible for 
the increase of the magnetic energy while the last term 
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is responsible for the destruction of magnetic energy due 
to diffusion. Since the equation for the magnetic field is 
linear it expected that after some transient behavior the 
amplitude of the magnetic field will grow or decay at an 
exponential rate 



b ~ b(x, t) e 



.It 



(2) 



where 7 is the growth rate and b is a bounded function 
in time. For steady velocity fields that will be examined 
here b is either time independent or a periodic function 
of time. 

The growth rate 7 and its dependence on the flow pa- 
rameters is the primary interest in this work. Given the 
functional shape of u the only control parameter in the 
system is the magnetic Reynolds number i?^ that in this 
work it is defined as 



R, 



U 



(3) 



r] is the magnetic difFusivity. U is the amplitude of the 
velocity field that is defined as 



U 



,1/2 



(4) 



where the angular brackets stand for spatial average. fc„ 
is the velocity inverse length-scale that we define through 
the vorticity of the flow w = V x u as 



,1/2 



/U. 



(5) 



Exploring the dependence of 7 on has been the 
subject of extensive research. For sufficiently small R^^^ 
the diffusion term in Eq.(IT]) will dominate and mag- 
netic energy will decrease exponentially with decay rate 
—7 ^ UkuR~^ (for <C 1). As i?^^ is increased the 
stretching term becomes important and above a critical 
value a flow can become an effective dynamo (7 > 0). 
This value of R^^ will be referred to as the critical mag- 
netic Reynolds number and will be denoted as R^ic- 
Finding the flow that minimizes Rj^j^^ is of importance for 
laboratory dynamo experiments on account of R^^ being 
an increasing function of power consumption which is an 
increasing function of cost. 

For large values of the magnetic Reynolds number the 
problem becomes increasingly complex with the number 
of degrees of freedom involved increasing like i?3/2. Due 

o ^ A/ 

to this complexity there is no general analytic way to 
estimate the growth rate of a dynamo (with the excep- 
tion of some special cases). Nonetheless, anti-dynamo 
theorems [T^, [l^ developed in the last century and up- 



per bounds on the growth rate |17H20l | have been proven 
useful in excluding certain classes of flows from giving 
dynamo action or restricting the scaling of 7 with R^^^ . 
From the Anti-dynamo theorems two important results 
that are relevant in this work are mentioned here. 

First, flows with only two non-zero components of the 
velocity field can not result in dynamo for any value 
of Rj^j [1^. Thus, for these flows there is no critical 



Reynolds number. In the present work we will refer to 
these flows as 2D flows (even if there is spatial depen- 
dence in the third direction). 

Second, time independent flows for which all three 
components of the velocity {ux,Uy,Uz) are non-zero but 
depend only on two of the spatial components (say x, y) 
can result in dynamo, but due to the absence of chaotic 
flow lines the dynamo growth rate will tend to a non- 
positive value as Rj^j tends to infinity [2l|, . For these 
flows it is thus expected that 



7 = o(f) Uku, 



(6) 



where the symbol "0(f)" stands for "smaller than order 
one" . This dependence however can be a very slowly 
decreasing function of i?^ [2^. These flows are referred 
to as "2,^0 flows, and the resulting dynamo is referred to 
as a slow dynamo. 

However, besides these classes of flows typical three- 
dimensional flows with a complex streamline top olog y are 
expected to be dynamos at infinitely large R,^j 2J]. For 



such flows the growth rate will approach a value that de- 
pends only on the amplitude, length-scale and structure 
of the velocity field and not on the magnetic diffusion rj 
ie : 



7 = 0{l) Uk^. 



(7) 



where the symbol "0(1)" stands for "same order as one". 
Such flows for which the dynamo growth rate tends to a 
positive value as R^^, tends to infinity will be called fast 
dynamo flows. 

With these restrictions in mind a definition of "an op- 
timal flow" can be given. The choice of optimization will 
of course depends on the application in mind. For ex- 
ample, an optimal flow can be based on R^^jf^ or on 7 
leading to different answers. Here we will restrict to the 
following questions: Given a family of flows of fixed ve- 
locity amplitude U and length-scale fc„ (i) which member 
has the smallest critical magnetic Reynolds number i?^^ , 
(ii) given R^^ which member has the largest growth rate 
j/{Uku), and finally (iii) which flow leads to the largest 
growth rate in the limit i?^,^ — > 00. As it is shown later, 
these questions do not have the same answer. 

Finally we need to restrict the family of flows that 
are going to be investigated. Since the estimate of the 
growth rate 7 and Rj^,c needs to be done numerically 
addressing the questions above for a large family of flows 
is formidable even for present day computing. It is then 
preferable to restrict to smaller families that are however 
good candidates for fast dynamo action based on their 
properties. 

For a fast dynamo the role of chaoticity and helicity of 
the flow has been emphasized as important ingredients. 
Chaoticity, the exponential stretching of fluid elements, 
is a necessary ingredient a for fast dynamo |2ll, . How- 
ever, it is not sufficient. Time dependent 2D flows can re- 
sult in chaos (positive Lyaponov exponent) but can be ex- 
cluded from dynamo flows based on the first anti-dynamo 
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theorem mentioned here. The reason for this behavior is 
that the flow is not only required to exponentially stretch 
the magnetic field lines but it needs to also arrange them 
in a constructive way so that when they are brought ar- 
bitrarily close they can survive the effect of diffusion. A 
quantitative measure of these effects can be obtained by 
multiplying Eq.(IT]) by b, space averaging and dividing by 
(b'^) and finally using Eq.([2]) to obtain 



7 + 9t In 



(b ■ (Vu)b) 



7?- 



(b2) 



(8) 



Performing a time average and using the fact that b is 
bounded the last equation can be written as 



7 = 7s - 7d 



where 7s is the time average of the first term in the right 
hand side of Eq.® and expresses the injection rate of 
energy by stretching, while is the time average of the 
second term and expresses the dissipation rate. A con- 
structive flow then has large 7s and small 7^. This can 
be obtained by the stretch-twist-fold mechanism [25| that 
aligns the stretched magnetic field lines so that they have 
the same orientation. It is expected to be achieved most 
efficiently if the flow is helical. 

Helicity is the other ingredient that is expected to im- 
prove dynamo action. It is a measure of the lack of re- 
flection symmetry of the flow Q and is related to the 
linking number of the flow lines. Although in general it 
is not necessary [2^ , it has been thought to improve dy- 
namo action and it is required for dynamos [27h29| . 
It is also considered important for the the generation of 
the large scale magnetic flclds that are observed in the 
universe jsol - lssj . 

In this work we are going to restrict ourselves in a 
family of flows that is both fully helical and are known to 
have chaotic flow lines, namely the ABC flows. The ABC 
flows include a wide range of expected dynamo behaviors 
that covers 2D flows, slow and fast dynamos. Particular 
members of this family have been well studied for dynamo 
action and this allows for a comparison with previous 
results. This choice is rather restrictive since it is not 
known a priori if the optimal dynamo flow belongs in the 
family of ABC flows. However, they provide a tractable 
set of flows to examine and a good starting guess. 

The ABC flows are reviewed in detailed the next sec- 
tion. 



It is an incompressible periodic flow with 4 independent 
parameters A, B, C and fc„. The flow has the property: 



w = V X u = fc„u 



(11) 



for all values oi A, B,C,ku- As a result it is an exact 
solution of the Euler equations. It has been studied both 
for its properties as a solution of the Euler equation, its 
relation to chaos [33 - [36| and for dynamo instability |37l - 
|39| but only for limited values of the parameters. Here 
it is attempted to uncover the dynamo properties for the 
whole family. 

With no loss of generality we can restrict ourselves 
only to flows of flxed wavenumber ku and flxed velocity 
amplitude U = V A"^ + B"^ + . With this restriction 
(^9) and property [TT] the energy of the flow E = |C/^, the 



enstrophy of the flow 



and the helicity H of the flow 



1 



H = -(w • u) 

2^ ' 



2 



(12) 



(13) 



have a flxed value. 

For flxed kinetic energy the parameters A, B,C live on 
the surface of a sphere of radius U and can be parame- 
terized using the spherical coordinates -0, 0: 



A ;7cos(V'), 

B ~ U sin(7/;) cos(0), 

C — U sm{tp) sm{(j)) . 



(14) 
(15) 
(16) 



Using the symmetries of the ABC flow (see [35|) we can 
restrict the examined parameter space. The flow is in- 
variant under the transformations 

[A,z]^[-A,z-TT/k^], (17) 
[B,y]^[-B,y-n/k^], (18) 
and [C,x]^ [-C,x~ n/ku]. (19) 

These symmetries allow to restrict the investigation to 
only positive values of A,B and C and thus reduce the 
examined parameter space to the range [0,7r/2] for both 
angles (j) and V- Since there is no preferred direction 
between {x, y, z) the growth rate is also going to be in- 
dependent under permutations, e.g. {B,C) — s- {C,B). 
More precisely the flow is invariant under the transfor- 
mations 



III. THE ABC FAMILY 

The ABC flow is named after V. Arnold [I], E. Bel- 
trami and S. Childress [l3], and is explicitly given 
by: 



U:c ~ Asin{kuz) + Bcos{kuy), 
Uy = C sin(fctix) -|- Acos{kuz), 
Uz = B sin(fc„?/) -I- C cos(fc„x). 



(10) 
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FIG. 1. A sketch of the parameter space for the ABC family 
of flows. RF marks the location of the Roberts Flow. The 
dashed lines indicate the location where two of the three pa- 
rameters {A, B, C) are equal and at their intersection is the 
1:1:1 flow. The grey lines enclose the region where stagnation 
points exist. 



These symmetries allow to intcr-change the values of any 
of the three parameters A, i?, C and were used to improve 
the estimates of the growth rate from the numerical sim- 
ulations. Finally, the change k —k changes the sign 
of the hclicity of the flow. This change however does not 
alter the resulting growth rate of the dynamo. Thus only 
positive values of k are considered. 

Depending on the values of the parameters A, B, C the 
flow can have eight stagnation points where all three com- 
ponents of the velocity are zero. These points exist only 
if the square of each the parameters A, B, C is smaller 
than the sum of the square of the other two [s^ (ze their 
squares can form a triangle) . The importance of the exis- 
tence or absence of stagnation points was emphasized in 
(T4| . It was noticed that the developed magnetic struc- 
tures changed from from "cigar-shaped" in the presence 
of stagnation points to "ribbon-shaped" in their absence. 

Figure [T] demonstrates the examined parameter space 
in the spherical coordinates (V', 0). Some of the points in 
this graph represent flows of special significance that are 
described in what follows. For (-0 = 7r/2, = 7r/2), {ip = 
7r/2, = 0) and (-0 = 0) two of the three parameters 
A, B, C are zero (A = B = 0, yl = C = OandB = C = 
respectably). The flow corresponding to these points is 
2D and thus there is no dynamo, 7 < 0. 

For (0 = tt/2) we have A=0, for (0 = tt/2) B=0, and 
for (0 = 0) C=0; for these values of ("0, 0) for which 
one of the three parameters A, B, C is zero, the resulting 
flow is a 2iD flow and thus a slow dynamo. The Roberts 
flow is a special flow in this subset for which the two 
nonzero parameters are equal. It has been studied for 
slow dynamo action in [ll|, [l^, H^. It corresponds to 
the values = 7r/2,0 = 7r/4), (-0 = 7r/4, = 7r/2), 
(■0 = 7r/4, = 0) and in figure [T]it is marked as RF. 

Flows which have two of the three parameters equal 
have additional symmetries and as will be shown in 



FIG. 2. Color-scale plot of the finite time Lyaponov exponent 
\t in the ip, plane. Bright colors imply large values of the 
exponents, while black imply zero or close to zero values. The 
contour lines correspond to the levels \t= 0.02, 0.04, 0.06 and 
0.08. The time of integration was r = 2 • 10*. 
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FIG. 3. Plot of the finite time Lyaponov exponent At for 
= 7r/4. The time of integration was r = 10^. 



the result section they are important. These flows are 
located along the line = it/A for B=C, the line 
■0' = arctan(l/ cos(0)) for A=B and the line ip = 
arctan(l/ sin(0)) for A=C. These lines are shown by 
dashed lines in figure [1] and divide the space in six com- 
partments. Each of these compartments is equivalent to 
the others due to the symmetries in Eq. (|20l21l22p . Thus 
each of these compartments will have the same number 
of maxima and minima of the growth rate. 

When all three parameters are equal A=B=C the flow 
has the largest number of symmetries. This flow is the 
most studied one in the literature and it is going to be 
referred to as the 1:1:1 flow. It is obtained for (0 = 7r/4, 
-0 = arctan(-\/2)) and is located in the intersection of the 
dashed lines in the diagram. 

Finally the region of the parameter space for which the 
ABC flow has stagnation points is enclosed by the grey 
lines in figure [T] 

ABC flows arc known to be chaotic [Sj-lsa- Finite 
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time Lyaponov exponents provide a measure of chaos 
poj . The finite time Lyaponov exponent Ar(a;o) for a 
point Xq is defined as: 



Ar(a;o) 



1 



■In 



|<5x(r) 



|<5x(0) 



(23) 



where |(5x(t)| is the distance of two particles that at time 
T ~ they were placed infinitesimally close to xq. If the 
flow is not ergodic, not all initial points xq of a chaotic 
flow lead to A^- > 0. To measure thus Ar of the flow, 
an ensemble of initial points xo needs to be considered 
out of which only those that belong to the chaotic sub- 
set will lead to Ar > 0. Here, Xt{xo) was calculated 
for the ABC flows and for 8000 initial positions xq dis- 
tributed uniformly in the domain [0,27r)'^. The distribu- 
tion function of the measured Lyaponov exponents was 
constructed and Ar of the chaotic subset was determined 
as the location of the peek in the distribution function. 

In figure [2] a color-scale plot of Ar is shown for the 
{ijj, (j)) plane and in figure [3] the finite time Lyaponov 
exponents are shown for (p = 7r/4. It is worth noting 
that the Ar of the most symmetric flow 1:1:1 is a local 
minimum (see also [13), while the largest values of Ar 
appear for {(f> = 7r/4, ip ~ O.ISStt) and for (0 ~ 0.127r, 
■0 ~ O.IGtt) and the equivalent points by symmetry. 



IV. DYNAMO RESULTS 

The advection diffusion equation ^ was solved in a 
triple periodic domain of size L ~ 2Tr using a standard 
pseudo-spectral method and a third order in time Runge- 
Kuta [U \^ . The resolution used varied from 32'^ grid 
points for small values of Rj^j (< 20) up to 256'^ for the 
largest values, i^j,^ > 500. Each run was evolved for 
sufficiently long time until a clear exponential increase of 
the magnetic energy was observed and the growth rate 
was calculated by fitting. 

The last parameter that needs to be defined is the ratio 
of the box size L over which the magnetic field is allowed 
to evolve in, to the period of the velocity field 27r/fc„. Due 
to the periodicity the product k^L can only be integer 
multiples of 2tt. Here we are going to examine two cases 
kuL = 2tt where the two lengths are equal, and fc„L = An 
where the magnetic field can evolve on a larger scale. 



A. ABC, kuL = 2n 

First the case fc„L = 27r is presented. In figure 2] color- 
scale images of the measured growth rate are shown for 
six different values of i?^^ . Each figure corresponds to 200 
different dynamo simulations. Using 20 different values of 
■0 in the range [0,7r/2] and 10 for (j> in the range [0,7r/4]. 
The symmetries in Ea. (|20ll2^ were used to fill in the 
values of the growth rate on the whole domain and on 
a denser grid. In each panel bright colors correspond 



to larger growth rate. The thick white lines show the 
location of zero growth rate. The thin black lines indicate 
where the growth rate is 0.05J7fc„. The dashed black lines 
as in figure [T] show the location on which two of the three 
parameters A, B, C are equal. Finally, it is noted that 
the simulations in these runs were performed on 32'^ and 
grid-point meshes. 

As Rm is increased the first flows that result in positive 
growth rates are the ones with two of the three parame- 
ters equal, while the third is equal to zero. This can be 
seen in the top left panel of flgurc[3]for Rj^j — 10 where 
most of the parameter domain has negative growth rate 
except the small bright regions at the end of the dashed 
lines. These flows correspond to a Roberts flow and are 
slow dynamos as discussed in the previous section. Thus, 
although they are slow, at small i?„ they are the most 
efficient at producing a dynamo {ie the fastest). 

The next flows that become unstable are the flows for 
which two of the the three parameters arc equal but 
smaller than the third. Thus they lie on the dashed 
lines in the graph "opposite" the Roberts Flow. This 
can be seen in the right top panel of figure S] that shows 
the growth rate for R^^^ = 14.3. In terms of the angles 
they correspond to the values {ip ~ O.ITtt, (j> = n/A), and 
{tp ~ O.SStt, = 7r/4 ± 0.127r). The exact location of 
these new maxima is shifting slowly away from the cen- 
ter as the magnetic Reynolds number is increased. Note 
that this flow is very close to the flow for which the max- 
ima for the Lyaponov exponents in figures [2] and [3] were 
found. It is also close to the flow A = 5, B~C^2 that 
was investigated in detail in [l3| , for this reason this flow 
is going to be referred as the 5:2:2 flow. At this value of 
the Reynolds number the Roberts flow is still the fastest 
dynamo in the family. 

As the magnetic Reynolds number is increased further 
the most symmetric flow 1:1:1 also results in dynamo. 
This is shown in the middle left panel of figure H] for i?^^- = 
20. At this value of Rj^, the 1:1:1 is a local maximum, but 
with smaller growth rate than the 5:2:2 flow and smaller 
than the Roberts flow that is still the fastest. 

As Rm is increased further the 1:1:1 flow stops being 
a local maximum and transitions to a to a third-order 
saddle point (monkey saddle point). This can be seen 
in the middle right panel for which R^^, = 25. The local 
maximum of the 1:1:1 flow, that was present at Rm ~ 20, 
splits to three local maxima that move along the dashed 
lines away from the 1:1:1 case who's growth rate has 
decreased. The growth rate for the 5:2:2 flow and the 
Roberts flow continues to increase. 

For Rm = 33.3 (shown in the bottom left panel) the 
three local maxima that were initially located close to 
the 1:1:1 flow have moved sufficiently away that the 1:1:1 
flow stops being a dynamo. This corresponds to the no- 
dynamo window that was observed early on in (37j . 

After further increase of i?„ the 1:1:1 flow becomes a 
dynamo again (although not a local maximum anymore 
but still a saddle point). For R^^^ = 100 shown in the 
bottom right panel most of parameter space is resulting 
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FIG. 4. Color-scale images of the growth rate in the ■(/), (f) plane for kuL — 2tt and for six different magnetic Reynolds numbers. 
Bright colors indicate larger growth rate. The thick white lines are the contour lines zero growth rate. Thin black lines are the 
contour lines of growth rate 7 — 0.05Uku- The dashed lines indicate the as in figure [1] show the location for which to of the 
three parameters A, B, C are equal. 



in dynamo action, with only exception the small areas 
close to the 2D flows {tp = tt/2, cj) = 7r/2), {ip = 7r/2, 
(/) = 0) and (-0 = 0). The growth rate of the Roberts flow 
has started to decrease and the fastest dynamo is given 
by the 5:2:2 flow. At this value of the "topogra- 
phy" of growth rate in the parameter space has become 
much more complex, with new local maxima appearing 
between the Roberts flow and the 5:2:2 flow. 

For larger values of i?„ grids larger than 64^ arc 
needed and it is computationally too expensive cover the 



whole parameter domain. Instead the investigation will 
be limited to flows which lie along the dashed lines where 
most of the maxima are located. In figure [S] we show 
the growth rate for three different values of the magnetic 
Reynolds number with the smallest value being equal to 
the value used in the last panel of figure |4l The 5:2:2 
flow results in the fastest dynamo at the largest value 
of = 1000. At this value of R^^ the 5:2:2 peek 
has moved to -0 ~ O.IGtt. Furthermore, for Rj^,, > 500 
two new local maxima appear close to the 5:2:2 flow for 
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FIG. 5. Non-dimensional growth-rate as a function of -0 for 
4> = tt/A and for 3 different magnetic Reynolds numbers. The 
vertical dotted line marks the location of the 1:1:1 flow. 
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FIG. 6. The rates 7s (dashed line, triangles) and yd (solid 
line, diamonds) defined in Eq. ^ for the kuL = 2n case and 
for Rj^j — 1000. The vertical dotted line marks the location 
of the 1:1:1 flow. 



slightly smaller and slightly larger values of ip. The local 
maximum vifhich at small values of i?^,^ was located at 
the 1:1:1 flovi' appears to return close to the 1:1:1 point 
and thus the most symmetric flow comes close to a local 
maximum again. Finally the slow decrease of the growth 
rate of the Roberts flow can be observed. 

In figure[S]the two rates 7s (dashed line, triangles) and 
7d (solid line, diamonds) deflned in Eq.([3]) are shown for 
the largest examined i?„ = 1000. The difference between 
the two curves gives the growth rate. The ratio of the two 
curves shows the percentage of the injected energy that is 
dissipated. Thus a constructive flow (in the sense that it 
aligns magnetic field lines pointing in the same direction) 
is expected to have a small value of 7^ compared to 7^. 
The flows close to the 5:2:2 flow (tp = O.IGtt) that have 
the largest growth rates are more efficient not only due 
to the larger stretching rate 7s that does not vary a lot, 
but also due to the relatively small value of jd- In the 
range -0 = O.Itt to -0 = 0.257r half of the energy injected 
by stretching goes to magnetic field amplification. For 



values of tp out of this range, only a small fraction of the 
injected energy goes to field amplification while the rest 
is going to the small scales where it is dissipated. 

Beyond the = 7r/4 symmetry line other local maxima 
of the growth rate were detected although it was not 
feasible to cover the entire parameter space. Here it is 
just mentioned that a local maximum was observed at 
(1/; = 0.27r, (j) = O.Itt) with growth rate close to the 5:2:2 
flow 7 ~ 0.12. 



B. ABC kL = 47r 

The case kuL ~ Att is examined in this section. Despite 
the fact that the flow is the same as in the fc„L = 27r 
case, the results arc different due to the additional space 
in which the magnetic held is allowed to evolve. The 
extra space gives rise to new modes that can develop 
with different growth rates. Although the modes of the 
kuL = 27r case are still present and grow at the same rate, 
they are not necessarily the fastest. Considering that in 
a numerical simulation only the fastest growing mode is 
observed, in the fc„L = Air case the observed mode will 
then be at least as fast as the kuL = 27r case. 

As before, the first fiows that result in dynamo arc the 
slow dynamos of the Roberts flow for which two of the 
three parameters A, B, C are equal and the third is zero. 
This case is shown in the right panel of figure [7] for R^^^ — 
2.5. Note that in this case the dynamo instability appears 
at much smaller values of R^^j . It is also remarkable that 
the mode whose growth rate peaks for the Roberts flow, 
appears to continuously extend all the way to the 1:1:1 
flow that is a saddle point at this stage. 

As i?„ increases further, the 1:1:1 flow becomes a dy- 
namo whose growth rate is a local maximum in the {ip, (f) 
plane. This is shown in the top right panel of figure [7| 
that corresponds to Rm = 10. Note also that this is con- 
trary to the kuL = 2tt case for which the 5:2:2 fiow was 
the second flow to result in dynamo. In the k^L = Att 
case and for this value of i?^^- there is no observed local 
maximum close to the 5:2:2 flow. 

As i?„ is further increased the growth rate of the 1:1:1 
flow is increased. At R^j = 25 the 1:1:1 fiow exceeds the 
Roberts flow in growth rate and it is the fastest dynamo 
for all ABC flows. This can be seen in the bottom left 
panel of flgurc [T] This is somehow surprising since this 
flow was never the fastest in the kuL = 2tt case. 

At even larger R^^^ however the growth rate of the 1:1:1 
ceases to increase while the 5:2:2 becomes a local maxi- 
mum and obtains comparable values with the 1:1:1 flow. 
This is shown in the bottom right panel in figure [T] 

The growth rate for larger values of i?^^ was calculated 
only along the symmetry line (jj = tt/A. It is shown as a 
function of ?/; and for three different values of R^^ in figure 
[8l The smallest value of R,^, corresponds to the results 
of the bottom right panel of flgurc [71 As the magnetic 
Reynolds number is increased the growth rate of the 1:1:1 
flow is decreasing while at the same time the growth rate 
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FIG. 7. Color-scale images of the growth rate in the ip, (j) plane for kuL = 47r and for four different magnetic Reynolds numbers. 
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of the 5:2:2 flow is increasing. At the largest examined 
value of the fastest dynamo is given by the 5:2:2 flow 
with a growth rate ^/{kuU) = 0.16 which is larger than 
its growth rate in the A:„L = 27r case. 

As in the previous section we plot in figure IH] the two 
growth rates 7^ (dashed line, triangles) and 7^ (solid line, 
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FIG. 9. The rates 7s (dashed line, triangles) and yd (solid 
line, diamonds) defined in Eq. ([9| for the kuL = 47r case and 
for _Rjj- = 500. The vertical dotted line marks the location of 
the 1:1:1 flow. 



diamonds) for = 500. In this case the the ability of 
the 5:2:2 flow to align field lines reducing dissipation is 
even more pronounced. Only one fifth of the injected en- 
ergy is cascading to the dissipated scales while the rest 
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is going in the amplification of the magnetic field. Un- 
like the kuL = 2-K case the 1:1:1 flow is also being more 
"constructive" with less than half of the energy going to 
dissipation. 

It is also worth comparing the general behavior of the 
growth rate with the results in the previous section. Al- 
though the fastest dynamo flows appear at the same lo- 
cation their growth rates are different thus it is not the 
same dynamo modes that are observed in the two cases. 
Also, in the fc„i — An case the dependence of the growth 
rate on the flow at the large R^^j is less complex than in 
the kuL ~ 2tt case with less local maxima and a smoother 
in general behavior. These differences indicate that the 
box size plays an important role in the dynamo behavior. 



V. SUMMARY AND CONCLUSIONS 

In this work the entire family of ABC flows was exam- 
ined for dynamo action. The dynamo growth rate was 
calculated as a function of the magnetic Reynolds number 
and for two length scales kuL = 2tt and = 47r. 
The questions that this work was attempting to answer 
were: (i) which flow has the smallest critical magnetic 
Reynolds number i?^^ j (^i) given R^ which flow has the 
largest growth rate 7/ (Uku), and (iii) which flow leads to 
the largest growth rate in the limit R^j — )■ 00. Although 
these questions can be posed for a larger family of flows 
the ABC-flows constitute a first step in obtaining some 
understanding. 

For this restrictive perhaps family of flows the answer 
to the first question is a simple one: The Roberts flow 
results in dynamo for the smallest value of This is 
true for both the k^L — 27r and the fc„L = Air case. Thus 
in small Reynolds numbers a well organized flow can do 
much better than a rapidly stretching (chaotic) flow. 

At larger Reynolds numbers new dynamo modes be- 
came unstable and a number of bifurcations are observed 
that lead to a complex "topography" of the growth rate. 
As i?J^J was increased this complexity is further increased 
and more local maxima appeared. This is particularly 
true for the fc„L = 27r case, while for the kuL = An case 
a smoother behavior was observed. 

Inspecting the growth-rate for a large number of flows 
as is done in this work also gives a wider perspective on 
the dependence of magnetic eigenmodes of the flow on 
R^. Some of the observed dynamo modes of a given 
flow can be be related (by continuous transform) to the 
modes of different flows. For example for small values of 
the slowest decaying mode of the 1:1:1 flow is related 
to the dynamo mode of the Roberts flow (see figure [7] top 
left panel). Thus the various bifurcations that can be 
observed by looking the growth rate of a single flow, can 
be interpreted as shifting or enlargement of local maxima 
in this wider point of view. The no-dynamo window of 
the 1:1:1 flow is such an example, which is the outcome 
of splitting and shifting of the initial maximum at the 
1:1:1 point. 



Finally, for relative large i?„ the 5:2:2 flow (■0 ~ O.IStt, 
(/> = tt/A) has the fastest growing mode (from the ex- 
amined flows) in both cases {k^L = 2tt and = Air). 
However, if this continues to be true for even larger values 
of the magnetic Reynolds number cannot be concluded 
from the present data. R^^ = 1000 is still far from the 
Rj^, — > 00 limit as can be seen from the finite value of the 
growth rate of the Roberts fiow which is a slow dynamo. 
Furthermore, as noted at the end of section (jIV A[) . fiows 
that were not on the (j) ~ n/A symmetry line were found 
with growth rates similar to the 5:2:2 fiow. The increased 
complexity of the growth rate as R^ increases makes it 
harder to estimate the fastest dynamo flow. If this con- 
tinues, then the location of the fastest flow in the {ijj, cj)) 
plane might not converge to a single point in the limit 
Rj^, — > 00 and question (iii) might not even have an an- 
swer. 

On the other hand the Lyaponov exponents, whose 
value does not depend on R^^^ do show some clear max- 
ima, which gives hope that a fastest dynamo flow in the 
— > 00 limit exists. However, although a correlation 
of the growth rate with the Lyaponov exponents is ob- 
served, it is definitely not sufficient to explain the de- 
pendence of the observed growth rates, at least not at 
the examined values of R^^, . In particular it is observed 
that the flow with the largest growth rate is close to the 
flow with the largest Lyaponov exponent. Nevertheless, 
the general dependence of the growth rate and of the 
Lyaponov exponent on the flow is quite different, with 
local maxima appearing at different locations. 

In addition it was found that the 5:2:2 flow that lead to 
the fastest growing mode besides having large stretching 
rate it was also very efficient at organizing the magnetic 
field lines as to minimize the magnetic energy dissipation. 
Furthermore at the examined values of R^^ the kuL = 2tt 
and the k^L = Air cases showed significant differences al- 
though the magnetic field lines are advected by the same 
flows. Thus the growth rate can not be determined by 
the stretching statistics of the flow alone. If these differ- 
ences cease to exist at larger R^^^ is a question of future 
work. 

Besides investigating larger R^^ there many other obvi- 
ous extensions of this work. First it would be interesting 
to extend these results in a larger family of flows, that 
also include non-helical flows. Harmonic velocity fields 
could be such a generalization. 

A differently oriented approach would consider turbu- 
lent dynamos. In this case, instead of prescribing the 
flow, a body force would be prescribed and the flow would 
be allowed to evolve dynamically. In such a study differ- 
ent limits of the kinetic Reynolds number Re would lead 
to different results. In the limit 1 <C Re <C dynamo 
growth rates depend on the small velocity scales and are 
possibly universal. In the other limit 1 <C Rj^j ^ Re 
it has been shown that the large scale flow plays an im- 
portant role especially for Rj^^ near its threshold value 

Finally the properties of dynamos beyond the linear 
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regime, where our understanding is much more hmited, 
is also a problem of considerable interest. At the non- 
linear stage, both the saturation levels of the magnetic 
energy and the involved length scales (large or small scale 
dynamo) depend strongly on the large scale properties of 
the flow. Thus a systematic study of a large number of 
flows can be helpful in that respect. 

These issues are going to be pursued in the authors 



future work. 
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